function x_next = ptm(x, a, b, c)
% 一维多项式三角映射 (Polynomial Trigonometric Map)
%
%   x_{n+1} = a * ( x_n + b*sin(3*x_n) + c )^2
%
% 输入:
%   x       - 当前状态
%   a, b, c - 系统参数
%
% 输出:
%   x_next  - 下一时刻的状态

    x_next = a * ( x + b*sin(3*x) + c )^2;
end
